Search for: All records

Using a graph representation of RNA structures, we have studied the ensembles of secondary and tertiary graphs of two sets of RNA with Monte Carlo simulations. The first consisted of 91 target ribozyme and riboswitch sequences of moderate lengths (<150 nt) having a variety of secondary, H-type pseudoknots and kissing loop interactions. The second set consisted of 71 more diverse sequences across many RNA families. Using a simple empirical energy model for tertiary interactions and only sequence information for each target as input, the simulations examined how tertiary interactions impact the statistical mechanics of the fold ensembles. The results show that the graphs proliferate enormously when tertiary interactions are possible, producing an entropic driving force for the ensemble to access folds having tertiary structures even though they are overall energetically unfavorable in the energy model. For each of the targets in the two test sets, we assessed the quality of the model and the simulations by examining how well the simulated structures were able to predict the native fold, and compared the results to fold predictions from ViennaRNA. Our model generated good or excellent predictions in a large majority of the targets. Overall, this method was able to produce predictions of comparable quality to Vienna, but it outperformed Vienna for structures with H-type pseudoknots. The results suggest that while tertiary interactions are predicated on real-space contacts, their impacts on the folded structure of RNA can be captured by graph space information for sequences of moderate lengths, using a simple tertiary energy model for the loops, the base pairs, and base stacks. 

Base pairing complementarity is central to DNA function. G.C and A.T pair specificity is thought to originate from the different number of hydrogen bonds the pairs make. Quantifying how many hydrogen bonds exist can be difficult because water molecules in the surrounding can make up for or disrupt direct hydrogen bonds, and the hydration structures around A.T and G.C pairs on duplex DNA are distinct. Large-scale computer simulations have been used here to create a detailed map for the hydration structure on A.T and G.C base pairs in water. The contributions of specific hydration waters to the free energy of each of the hydrogen bonds in the A.T and G.C pairs were computed. Using the equilibrium fractions of hydrated versus unhydrated states from the hydration profiles, the impact of specific bound waters on each hydrogen bond can be uniquely quantified using a thermodynamic construction. The findings suggest that hydration water in the minor groove of an A.T pair can provide up to about 2 kcal/mol of free energy advantage, effectively making up for the missing third hydrogen bond in the A.T pair compared to G.C, rendering the intrinsic thermodynamic stability of the A.T pair almost synonymous with G.C. 

Performing full-resolution atomistic simulations of nucleic acid folding has remained a challenge for biomolecular modeling. Understanding how nucleic acids fold and how they transition between different folded structures as they unfold and refold has important implications for biology. This paper reports a theoretical model and computer simulation of the ab initio folding of DNA inverted repeat sequences. The formulation is based on an all-atom conformational model of the sugar-phosphate backbone via chain closure, and it incorporates three major molecular-level driving forces—base stacking, counterion-induced backbone self-interactions, and base pairing—via separate analytical theories designed to capture and reproduce the effects of the solvent without requiring explicit water and ions in the simulation. To accelerate computational throughput, a mixed numerical/analytical algorithm for the calculation of the backbone conformational volume is incorporated into the Monte Carlo simulation, and special stochastic sampling techniques were employed to achieve the computational efficiency needed to fold nucleic acids from scratch. This paper describes implementation details, benchmark results, and the advantages and technical challenges with this approach.